ON THE CHARACTERISTIC POLYNOMIAL 
OF THE ALMOST MATHIEU OPERATOR 
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Abstract. Let Ag be the rotation C*-algebra for angle 9. For d = p/q with p 
and q relatively prime, Ag is the sub-C*-algebra of Mq(C{T'^)) generated by a 
pair of unitaries u and v satisfying uv = e^^^^vu. Let hg x = u + u~^ + X/2(v + 
v~^) be the almost Mathieu operator. By proving an identity of rational 
functions we show that for q even, the constant term in the characteristic 
polynomial of hg^x is (-1)9/2(1 + (A/2)'?) - (zf + z^" + (X/2)i{z!^ + z^")). 



1. Introduction 



Let 9, A, and tp be real numbers with A positive. The second order difference 
operator Hg^\^^ on ^^(Z) given by 

H0^x,^{C){n) = ^(n + 1) + - 1) + Acos(27rn6' + V')^(ri) 

for ^ G £^(Z) is cahed the almost Mathieu operator. Hg,\^^ is a discrete Schrodinger 
operator which models an electron moving in a crystal lattice in a plane perpen- 
dicular to a magnetic field. 

An object of much study has been the spectrum ct(0, A) — U^a{Hg^x^^). In 
Hofstadter calculated a{9, 2) for 6 — p/q and 1 < p < q < 50. The remarkable 
pattern he found is called Hofstadter's butterfly. For irrational 9, a long standing 
concern has been the connectedness and Lebesgue measure of a{9, A) and the la- 
belling of the gaps, about which quite a bit is now known (see jAjj . [Tk. . and 
for spectacular recent advances as well as jAvMsj . [55] . jB], jcEvj . [T7r] for earlier 
work). In addition there has been numerical work on computing the spectrum to 
high accuracy for large q |Ai| A2 , l] . 

Let Ag be the rotation C*-algebra (see (b)). For 9 ~ p/q with p and q relatively 
prime and p = e^'^^ let 



Ug 



1 

0/ 



and vg — 



(P 



V 



V 



i.e Ug cyclically permutes the elements of the standard basis and vg is a diagonal 
operator. Then define u,w : ^ Mg(C) by u(zi,Z2) = ziug and v{zi,Z2) = Z2Vg. 
Then uv = pvu and Ag is the C*-algebra generated by u and v (see [B]). The 
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operator he^x = u 



X/2{v 



contains all the spectral information of 



He,\,^ in that Sp{hg^\) = ag^x := LI^Sp{Hg^\,jp). 

The main tool in the analysis of ag^\ is A.g^\, the discrete analogue of the dis- 
criminant. For = p/q, Ae,\{x) = Tt{Ai{x) ■ ■ ■ Aq{x)) where 



Ak{x) 



X — Xcos{2TTkp/q + n/ {2q)) — 1 
1 



Below are the first few values of this polynomial. Note the form of Ag^x so 
displayed depends only on the denominator q; however, ^g = 2 cos(27rp/(7) depends 
on the numerator p. 

Ae,2{x) for 9 = p/q and = 2cos(27re') 



2 




4 


3 


x^- 


6x 


4 


x"- 


8x^+4 


5 


x'' - 


lOx^ + 5(3 - 


6 


x""- 


12x'' + 6(5 - 


7 


x' - 


14x5 _ 


8 


x^~ 


16x^ + 8(9 - 


9 


x'- 


ISx'^ + 9(11 


One can 


calculate fo 



conjectural (from numerical evidence). A deeper understanding of the structure of 
Ag,\ would be quite interesting. 



k 


coefficient of x'' '^'^ in Ag^\ {fj, = 


A/2) 


1 


-g(i + M') 




2 






3 


~<lil^s{''f)f^' + + {"-/)- i<l - 






+ (1+(V) -(9-6)?e + 6fl)M' 





The connection with the characteristic polynomial of hg^x is given by 

(1) det(x/, - heMzi,Z2)) = Ag^xix) + zf + zp + (A/2)«(z| + z''^) 

and thus (T8,A = Ag\[-2{l + {X/2)i), 2(l + (A/2)9)]. Indeed, Ag^x{x) can be written 
as a determinant (c.f. Toda T, §4]) 



/ai 
1 



(2) 



^p/q,\i^) = det 



Vi 



1 

a2 



1\ 



+ 2 (-l)'' + (A/2)« 



where all the other entries are and ak — x ~ Acos(27rfcp/g + ■n/{2q)). Since 

(3) = {-ir\/,M^) 
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the coefficient of ai^^f^fe+i) jg Q for < /c < q/2. 

The main result of the paper asserts that for ai = 2 cos{2ttIp / q) and 1 < fc < q/2 
we have 







where the summation is over ah subsets of {1, 2, 3, . . . , g} obtained by deleting 
k pairs of adjacent elements - counting 1 and k as adjacent. This is proved by 
establishing the following identity for fc > 3 and g > 2fc — 1 



g-2(fc-l) 



J. J. Cj.— ij — 1 _|_ +1) 

41 = 1 ifc=ifc-i+2 j = l ^ ' 



2k-2 
i=fe+l 



fe-2 

)n( 

1=1 i=-l 

g-2 



11=3 Jfc_2=U--3+2 j = l 



(a; ^ + x^) ^{x 'I + x'i) X - X - 1 r \;x --^x• ■'j - 



We then use this to show that for a; = 2 cos{2Trlp / q) 
(ai I 1 \ 



det 



1 a2 



Vi 



■■. 1 

1 aj 



g = (mod 4) 
= {A g = 1,3 (mod 4) 
8 q = 2 (mod 4) 



From this we show that the constant term (i.e. the coefficient of x^) in det(a;/g — 
he^x{zi,Z2)) is 

(-1)^/^2(1 + (A/2)^)) - (z? + zr" + {m'{4 + ^2')) 

when q is even. When g is odd it follows from Q that the coefficient of a;" is 

-(z? + zr + (A/2)n4 + ^2"'))- 

Similar, though simpler, reasoning shows that the coefficient of x'' ^ is —q{\ + 

A/2) and that the coefficient of x"?-* is (A/2)4g(g-3)/2+(A/2)2g(g-4-2cos(27r6l)) + 

q{q-i)/2. 



2. The Main Theorem 



Let us use the following notation: let ai , . . . , a„ be elements of a commutative 
ring and let 

ai 1 
1 a2 1 



(ai, 02, . . . , a„) 



as 



1 a„ 
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and 



lai, a2, . . . , a„ 



ai 
1 



1 

02 



0-3 



1 



1 1 an 

The first matrix is a tridiagonal matrix with I's on the sub and super-diagonal and 
O's elsewhere. The second matrix is the same tridiagonal matrix with in addition 
I's in the upper right and lower left corners, all other entries are 0. Expanding 
along the bottom row we have 



(4) [ai,a2,...,a„] = (ai,a2,...,a„) - (a2,a3,...,a„_i) + 2(-l)" ^ 
and 

(5) [-ai, -02, . . . , -an] = (-l)"[ai, a2, • ■ ■ , a„] + 2(-l)"-i. 
Rewriting equation (j^J we have 

(6) Ap/,a(x) = [ai, . . . , aj + 2((-l)« + (A/2)') 

Notation 2.1. (i) For < fc < n/2, let S[l] = {/ C {l,2,...,n} | |/| = 

n — 2k and / is obtained from {1, 2, . . . , n} by deleting k pairs of adjacent 
elements}. S['^] = {0}, ^ {{!}, {3}, {5}, . . . , {2fc + 1}} , 

S[;]= {{1,2,3,..., n}]. 

(ii) For < /c < (n - 1) /2, let S' [I] = {/ C {2, 3, . . . , n} | |/| = n - 2fc - 1 and 
/ is obtained from {2,3,..., n} by deleting k pairs of adjacent elements}, 
yp^+i] = 0, S'[''+'] = {{2},{4},{6},...,{2fc + 2}}, S'[-] = 
{2,3,. ..,n}. 

/zmJ For S a collection of subsets of {1, 2, . . . , n — 1} let S V {n} = {/ U {n} \ 

leS}. 

(iv) For < A: < n/2, let S'[g = {/ C {1, 2, 3, . . . , n} | |/1 =n-2k and / is ob- 
tained from {1, 2, . . . , n} by deleting k pairs of adjacent elements, counting 
{n,l} as an adjacent pair}. S[^^] = 0, ^P\t'] = {{1}, {2}, . . . , {n}}, 
^K] ={1,2,3,. ..,n}. 
(v) For 01,02,03, . . . , a„ elements of a commutative ring, and / = {11,^2, *3, 
. . . , ifc} C {1, 2, 3, . . . , n}, let a/ = ai^^ai^ai^ ■ ■ ■ ai^ . We shall adopt the 
convention that 00 = 1 

Part (ii) of the next proposition goes back to Sylvester's original paper on con- 
tinuants IS); part {iv) is a straightforward extension of this. For the reader's con- 
venience we present a proof. 

Proposition 2.2. (i) Suppose 1 < k < n/2, then S[l\ = (^["j:^] V {n}) U 

[n/2] [(«-l)/2] [(«-2)/2] 

fc=o j(.s[l] k=o ies["-'] k=o ies["-^] 
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(iv) When n is odd, 

[n/2] [n/2] [(«-l)/2] 

E(-i)'= E «^ = E(-i)' E E (-1)'= E «^ 

fc=o fc=o /gsH fe=o /6S'["-] 

When n is even, 

[n/2] [n/2] l{n-l)/2] 

E(-i)' E «^ = E(-i)' E ^^^ + (-1)"^'- E (-1)' E «^ 

Proo/. ("i; Let / e • If n ^ / then n - 1 ^ / and so / e S'[feZi] • Suppose n G J. 
Let K = {1,2,3,..., n} \/ and / = / \ {n}. Then / = {1,2,3, .. .n - 1} \ K; so 
i e . Hence / = / U {n} 6 V {n}. 

('zzj Let us assume that n = 2m is even. The same idea works for odd n but the 
proof is shghtly simpler. Observe 

[n/2] m-1 

E(-i)' E '^^^ E '^/ + E(-i)' E «/ + (-ir 

m — 1 771—1 

= ^ a, + a„^(-lf 5] + (-1)- 

^ m — 1 m—1 

= «" E «^ + E(-i)'E «^ +E(-i)'E «/ + (-ir 

m—1 m — 1 

= «nE(-i)'E «^-E(-i)'E «^ 

[(n-l)/2] [(«-2)/2] 

= E (-i)'E E (-i)'E 

I'm; For I e let ii'i = {1, 2, 3, . . . , n} \ / and K2 = {2,3, . . . ,n}\I. min{i | 
i e Ki} is odd if and only if 7 G S'[^] and min{i | i G K2} is even if and only if 

(iv) Suppose n = 2m. Then 

[n/2] m-l 

^(-1)'= E E «^ + E(-i)' E «/ + (-ir E «^ 

Cm \ m 

m m — 1 

= E(-i)' E «/-E(-i)' E «/ + (-ir- 

fe=o 7es[^] fe=o /es'f";^!] 
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The case of n odd is similar. 



□ 



Corollary 2.3. Let ai, 02, . . . , a„ be elements of a commutative ring. 

[n/2] 

(i) (ai,...,a„) - ^(-1)'^ J2 



(n) 



[n/2] 
[„/2] 



odd 



Proof, (i) For n = 1 the left hand side and the right hand side equal ai. Both sides 
satisfy the same recurrence relation. 
(ii) By equation JHJ 

[ai, . . . , a„] = (ai, . . . , a„) - (02, . . . , a„_i) - (-1)"2 

so the result now follows from (i) and Proposition 12. 21 {Hi). □ 

Proposition 2.4. Let 1 < p < q be relatively prime, 9 = p/q, and ak — '2 cos(27rfc6'). 

{0 q = (mod 4) 
2 g = 1, 3 (mod 4) 
-4 (7 = 2 (mod 4). 

Proof. Let be the qth Chebyshev polynomial of the first kind. The constant 
term of Tg is for q odd and (—1)'/^ for q even. The result now follows from the 
identity (see e.g. |E] §1.2]) 

Y[{x-a,) = 2iTq{x/2)-l). 

□ 

The statement of the main theorem follows. Its proof will be given at the end of 
the next section. 

Theorem 2.5. Let I < p < q be relatively prime, Ok = 2cos(27rfc0), and 9 — p/q. 
Forl<k< q/2, 

^ aj = 0. 

Corollary 2.6. Let \ < p < q be relatively prime, 9 = p/q, A > 0, and Ok = 
Xcos{2TTk9). Then 

{0 q = (mod 4) 

2(1 + (A/2)'') (7=1,3 (mod 4) 
-4(1 + (A/2)«) g EE 2 (mod 4) 

and Ag^x{0) = (-1)9/22(1 + (A/2)9) for q even. 
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Proof. Suppose q is even. By Theorem 12. 51 all the terms of 

[9/2] 

k=o ies[l] 

are zero except the terms for fc = and k — q/2. The term for fc = is aia2 ■ ■ • a^. 
The term for k — q/2 is (—1)'/^. Thus when q — 4m we have by Proposition [^| 

[ai, 02, . . . , a,] - aia2 • • • a, - (-1)''2 + (-l)«/22 = 0, 

and when q = 4m + 2, 

[ai, 02, . . . , aj = aia2 • • • a, - (-1)«2 + (-1)^/22 = -4(1 + (A/2)«). 

To obtain the final claim we apply equation ®. □ 
From the corollary and equation we have the theorem which corrects an error 
in I^eTI p. 232]. 

Theorem 2.7. The coefficient of in det(a;/q — hg \{zi., Z2)) is 

- {zl + + (A/2)nz| + z-'')) + (-1)^/22(1 + (A/2)«)) 
w/ien q is even and ~{z\ + z^'^ + (A/2)'(z| + z^ ^)) when q is odd. 

3. Proof of the Main Theorem 

Theorem 3.1. Suppose 01,02, . . . , o^ are elements in a commutative ring and let 
Og+i = oi. For / C {1,2,..., g}, let I"^ — {1,2, . . . ,q} \ I, be the complement of I 
in {1, 2, . . . ,q}. Then 

q-2(k-l) q-2(k~2) q k 

E E E ■■■ E 11"*. "^.+1 

-g-2(fe-2) 9-2(/c-3) g-2 fc-2 

-0102 ^ E "' E n°*j"*j + l "'/"'J+l- 

ii=3 i2=il+2 ife-2=»fc-3+2 j=l 

Proof. Recall that elements of 5 [^] are obtained by deleting k adjacent pairs {i,i + 
1} from {1, 2, . . . , g}, counting q and 1 as adjacent. So if G S'['] then / = 
{«i,ii,«2, j2, ■ • • ,ik,jk} with 1 < ii,ji = ii + 1 < «2, ■ ■ • , jfc-i = «fe-i + 1 < ifc < g 
and either jk = ifc + 1 if «fe < 9 or jk ~ 1 if i^ = q. 

Now let T[^] = { {ii,ii,i2,j2, • ■ ■,ik,.ik} \ 1 < ii,ji = ii + 1 < 12, • ■ -Jk-i = 
ife_i + 1 < ifc < g, jfe = ifc + 1 }. Define (j) : {1, 2, . . . , g, g + 1} {l,2,...,gf} by 
(f>{q + 1) = 1 and (f>{i) — i for i < q. Then 0^(7) = o/ for ^ G . 

If / = {ii, ji,«2,j2, • • ■ ,ik,jk} and ii = 1 and ik ^ q then ^ S[l] because 

4>{jk) = (/"(ii) — 1 and the pairs must be disjoint. So let = { {1, 2, ii, ji, . . . , 

ik-i,jk-i,q,q+ 1} I 3 < ii, ji = ii^+ 1 <i2,.. .,ik-i <q - 2,ifc-i = ik-i + 1 }■ 

For / e T[l] \r[l\, e and : Tg] \T'[a ^ ^[^] is a bijection. 

This with the identity 0^(7) = 07 proves the theorem. □ 

Lemma 3.2. (^ij For q>l, 

2J,X +X) [X +X ) ~ (^^-2_2.2)(^-g-l+^,+ l)- 
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(ii) For k >\ 



i—1 ^ ^ ^ ' 

Proof, (i) One checks directly that the formula holds when q = I, then (i) follows 
by induction on q. 

(ii) follows from the identity 

— 1 



— x'^^){x~^~^ + x^+^) {x~^ + x^){x~ 



□ 



Corollary 3.3. For q > 5 

(x-1 + x'i)(x-''-^ + Xl+^ ' ^ ' 

{x-^ ~ x^){x-i^'^ - x'l-^) 

~ _ a;4)(a;-6 _ x6)(a;-'?+i + + x'i){x-i-^ + i-^+i) 

Proof. By Lemma 13.21 (i) 

+ -^ + x')-^{x-'-^ + x'+^)-^ 

i=3 



_ x2-^(^jr-q+l _|_ ^.g-l) (3,-2 _ x'^){x-^ + X^) 



-3 



The result then follows by multiplying both sides by (x ^ + x){x ^ + x'^){x 
Theorem 3.4. For k > 1 and q > 2k — 1, 

g-2(fe-l) g-2(fe-2) q k 



(7) 



Proof. We prove the equation by induction on k. When fc = 1 the equation holds 
by Lemma 13.21 (i). Lemma 13.21 (ii) shows that for arbitrary k the formula holds for 
(J = 2fc — 1; so we fix k and proceed by induction on q. Let Sk.q and T^^q denote 
respectively the left hand and right hand sides of equation 

If we write Sk,q as a sum of two terms, the first in which ik < q and the second 
when ik — q, we see that Sk,q satisfies the recurrence relation 

Sk,q = Sk,q-i + + + x''+Y^Sk-l,q-2 
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Thus we have only to show that Tk^q satisfies the same relation. Now 



'''''' nlii^-'' - ^'') uto(^-^'-'^ + ^'-') 

and 



The proof of the recurrence relation for Tk^q is thus reduced to verifying that 

(^^-{q-{k-l) _ ^q-(k-l)-j^^-{q-(k-l) _^^q-(k-l)'^ 
^-{q-{2k-l)) _ ^q-{2k-l) ^-2k _ ^2k 



□ 



x-'i+xi {x-'i + xi){x-'J-^ + X<i+^) 

Theorem 3.5. For k > 3 and q > 2k — 1, 

{x-^ + xY\x-^ + x^)-\x''' + x'^y^x-^-^ + 

g-2(fe-2) (}-2(fe-3) q~2 k-2 

il=3 12=11+2 ifc-2=ifc-3+2 j = l 

^ - x'^-'){x-'' - x^) Utllli^-'^' ~ ^'-') 

Proof. Let us denote the left and right hand sides of the identity by 5"^^^ and T^^q 
respectively. By Corollarv l3.3l 5-; = T^^q. We write 5^^, as the sum of two terms: 
in the first ik-2 < q ~ 2 and in the second ik-2 — q — 2. As in the proof of the 
previous theorem we obtain a recurrence relation, in this case; 

Sk.q = Sk,q-i{x-''-^ + + x"-') 

+ Sk-i,q-2{x-' + x'')-\x-'i-^ + 

It is routine to verify that Tk,q satisfies the same recurrence relation. □ 
Subtracting equation (jHl from equation ((JJ) yields. 

Corollary 3.6. 

q-2{k-l) q-2{k-2) q fc 

E E ••• E \\(xr'^'+x^r'{xr''~'+x^'^'r' 

11=1 i2=ii+2 ifc=ii._i+2 j=l 

- + x^)-\x-^ + x^)-\x-'^ + + x'^+^)-^ 

q-2{k-2) g-2(fc-3) q-2 k-2 

X E E •■• E X{{x-'' +x'^)-\x-''-^ +x'^+^)-^ 

il=3 12=21+2 ife-2=»fc-3+2 j=l 
2fc-2 

i=fc+l 



(9) 



fe-2 



\{(x'^ - x-'^) n( 

k=l 
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Proof of theorem \2.5\ We recall that I < p < q and p and q are relatively prime. 
We set 9 — p/q and aj = 2 cos(27rj'6'). We shall split the proof into two cases. 

Case 1: q ^ (mod 4). When q ^ (mod 4) aj / for all j; moreover when 
X = e^'^^^ , x^"^ ^ 1 and x^'-''"*^ ^ —1 for all i. Thus the denominator on the right 
hand side of © does not vanish but the numerator does. Hence by Theorem 13. II 

Upon multiplying by aia2 ■ ■ ■ Oq we obtain that 

ai = aia2---aq ^ (a/c)"^ = 0. 
ies[l] ies[l] 

Case 2: q = (mod 4). Again we wish to show that X]/gs[''] — ^^^d so we must 
multiply both sides of equation @ by rii=i(^~' + evaluate at x = e^'^*^. 

The denominator of the right hand side of © is zero when x'^'^ = 1 or x^^'^^^'> = 
— 1, i.e. when i — j — q/A; the corresponding factors are x^"^/^ — x''/^ and x^'^''/'* + 
x^i/"^ respectively. 

Apart from the factor x^'' — x'', the numerator of the right hand side of equation 
® is zero only when x^*^'^*-* = 1, i.e. when i = q/2. This produces the factor 

X-g/2 _ 2.9/2 

which cancels one of the zeros in the denominator. The other zero is 
cancelled when we multiply by n?=i(^~*+^*)- Hence the product of YU^i {x^'^+x"^) 
and the right side of Q is zero when x — e^'^'^. □ 
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